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CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims priority from US Provisional Patent Application No. 60/395,715 
filed 07/12/2002, which is incorporated herein by reference. 

FIELD OF THE INVENTION 

This invention relates to computer-implemented methods within the field of finance. More 
specifically, the invention relates to methods directed at valuing complex uncertain cash flow 
streams and determining the best way to hedge exposure to such streams. 



BACKGROUND 

Frequently it is desirable to assign a price to an asset that is not currently traded, or to 
independently verify a price of an asset that is traded. There are two standard methods 
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for assigning such prices (both of which earned Nobel prizes). The first is the capital asset 
pricing model (CAPM) [1] which is a method for assigning a price to an asset that has a 
duration of one period of time (such as one year) and is not traded at intermediate times. 
The second is the Black-Scholes method [2] for pricing financial derivatives, which assigns a 
price to a derivative (such as an option) on a underlying asset that is traded. 

A derivative security is a security whose payoff is determined by the outcome of another 
underlying security. For example, a stock option on a traded stock is a derivative, since the 
final payoff of the option is completely determined by the value of the stock at the terminal 
time. The Black-Scholes equation is the standard method for determining the value of 
derivatives. It is based on the fact that, in continuous time, it is possible to replicate the 
payoff of a derivative by a portfolio consisting of the underlying security and a risk free asset 
(such as a U. S. Treasury bill) with a fixed risk free interest rate r. The fractions of the 
portfolio devoted to each of its components is adjusted continuously so that the portfolio's 
response to changes in the underlying security will perfectly mirror the response of the 
derivative. This adjustment process is termed a replication strategy. It is argued that the 
value of the derivative is equal to the cost of replication strategy. The cost V is determined 
by the Black-Scholes equation, which gives the cost (or value) V(x,t) for values of x > 0 
and 0 < t < T where T is the terminal time of the derivative. Specifically, the equation is 

rV{x, t) = V t (z, t) + V x {x, t)rx + \V xx {x, t)x 2 a 2 . (1) 

Here x denotes the value of the underlying security, r is the annual risk free interest rate, 
and a 2 is the annual volatility of the underlying security. The notations V t ,V x ,V xx denote, 
respectively, the first partial derivative with respect to £, and to x, and the second partial 
derivative with respect to x. The equation is solved with the boundary condition V(x,T) = 
F(x(T)), where F denotes the payoff at time T of the derivative. For example, if the 
derivative is a call option with strike price K, then F(x(T)) = max (x(T) - K, 0). 

A more general situation is where a payoff depends on a variable x e but this variable is 
not traded. For example, an option depending on a firm's revenue is of this form, because the 
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revenue (which is x e in this case) is not traded. Hence the payoff depends on a non-traded 
underlying variable. In these situations it is impossible to form a replicating strategy using 
the conventional Black-Scholes equation because it is impossible to trade the underlying 
variable. The assumption underlying the conventional Black-Scholes equation breaks down. 

Such problems have been studied by other researchers. The idea of using a market hedging 
strategy to minimize the expected squared error between the final value of the hedge and the 
actual payoff was proposed by Follmer and Sondermann [3] , who showed that it was possible 
in principle. Because minimizing the expected squared error is equivalent to orthogonal 
projection of the payoff onto the space of marketed payoffs (under a standard definition of 
projection), the method is often referred to as projection pricing. Their analysis, however, 
is purely abstract and does not exhibit any practical method for explicitly finding the hedge 
by solving a partial differential equation or a discrete version of it. 

Merton (in his Nobel Prize acceptance speech [6]) emphasized the importance of the 
problem. He proposed a procedure based on the original Black-Scholes equation, but it 
is essentially an ad hoc method that does not coincide with projection. He and Pearson 
[7] studied a general framework of incomplete markets and formulated prices based on a 
consumer maximization problem. Again their method is abstract and does not provide a 
direct formula for the value that could be used in practice. 

Some practitioners have proposed various ad hoc translations of the standard Black- 
Scholes equation to specific situations. For example, it is common practice to artificially 
increase the volatility of the underlying variable in an attempt to recognize that the variable 
is not really a traded asset. These methods are not based on optimality, do not fundamentally 
revise the original Black-Scholes equation, nor have any other real theoretical basis. 

An approach that directly addresses the problem presents a solution in terms of a market 
price of risk, which is applicable to all derivatives of a non-marketed variable. The market 
price is difficult to measure, but it has sometimes been estimated [8]. However, even if the 
market price of risk is known, it does not lead to a hedging strategy. 
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Overall, there has not been an effective and practical method proposed that prices, opti- 
mally hedges, and computes the residual risk (after hedging) of derivatives of non-marketed 
variables. 

SUMMARY OF THE INVENTION 

Embodiments of the invention provide computer-implemented methods, based on projection, 
for determining the value of a derivative of a non-traded variable. In addition, some em- 
bodiments provide a method for determining the optimal hedging strategy and its residual 
risk. 

Several embodiments of the invention are centered on a new, extended Black-Scholes 
equation of the form 

rV{x € ,t) = V t (x e ,t) + V Xe (x e) t)x e [ti e - /3 em (fi m - r)} + \V XeXe {x e ,t)x 2 e al, (2) 

where 

Am ~ Pem&e/ ®m ■ 

The underlying variable x e is not marketed, (or equivalently, not traded). By this we 
mean, throughout, that it is not marketed as a security that can be held without cost and 
pays no dividend while held. Thus, although oil can be bought and sold, it is not traded 
in this sense because there are storage costs. However, a futures contract on oil is traded 
as a security. In some cases, a variable that is not traded in our sense can be converted to 
a traded asset; (for example, a dividend-paying stock can be so converted by subtracting 
the present value of its dividend). Such cases do not require the methods of this invention 
(although the methods of this invention may lead more quickly to the correct result). If our 
methods are applied to a case where the underlying variable is in fact traded, the extended 
equation will reduce to the standard Black-Scholes equation. 
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The extended equation (2) uses, in addition to the underlying variable :r e , an asset x m 
that serves as a market representative. This representative is itself a marketed variable, or 
a linear combination of marketed assets, in the sense described above. This variable enters 
the equation indirectly through specification of certain parameters of the equation. The 
constants fx e and \i m are the drift rates of x e and x m , respectively. p em is the correlation 
coefficient between x e and x m . a\ and are the annual variances of x e and x m , respectively. 

The replacement of r by \fi e -PemiPm— r)] as a coefficient of V Xe (x e , t)x e in the partial dif- 
ferential equation distinguishes this extended equation from the conventional Black-Scholes 
equation. This coefficient accounts for the fact that perfect replication is not possible in all 
cases. 

The market representative x m can be chosen in several ways, each of which leads to 
an identical result. In some embodiments x m is chosen to be the Markowitz (or market) 
portfolio of risky assets. According to other embodiments, x m is chosen to be a portfolio of 
marketed securities most correlated (related) to the underlying variable x e . This method is 
simpler than choosing x m to be the Markowitz (market) portfolio and furthermore leads to 
an optimal mean-square hedging strategy. 

Embodiments of the invention include an implementation based on a solution to a 
continuous-time version of the extended Black-Scholes equation. Thus, according to an 
embodiment, a method is provided for pricing a financial derivative of a non-marketed vari- 
able represented by x e . The method includes determining a market representative x m that is 
useful in determining the value of the financial derivative. Information associated with the 
non-marketed variable and the market representative is retrieved. A solution to an equation 
involving a variable V(x e , t) representing a price of the financial derivative is then calculated. 
The equation, which in this embodiment is preferably the extended Black-Scholes equation, 
includes a coefficient involving the information associated with the market representative x m 
and the non-marketed variable x e . An output including a calculated price of the financial 
derivative is then generated. Another embodiment includes a solution using formulas based 
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on risk-neutral processes. 

The price function can also be found in a discrete-time implementation. The underlying 
process can be discretized directly, or a finite-state model can be used. For example, in an 
embodiment of the invention, a method is provided for pricing a financial derivative of a non- 
marketed variable represented by x e . The method includes retrieving information associated 
with a suitable market representative x m . A solution is then found to a system of equations 
involving a variable V(x ei k) representing a price of the financial derivative at discrete times 
indexed by k. The equations include a coefficient involving the information associated with 
the variable x m associated with the market representative. An output including a calculated 
price of the financial derivative is then generated. 

In another embodiment, a computer-implemented method is provided for pricing a fi- 
nancial derivative of a non-marketed variable represented by a finite-state variable B. The 
method includes determining a market representative, represented by a variable A. Risk- 
neutral probabilities are then calculated using a binomial lattice model associated with the 
non-marketed variable and the market representative. Values of a variable V on the lattice 
corresponding to the variable B are also calculated. The calculated values of V represent a 
calculated price of the financial derivative. 

Embodiments of the invention may also involve computing an optimal hedge and other 
variables of interest. The above embodiments may be easily implemented in various different 
ways on any type of computer, and may be realized as instructions stored on a computer- 
readable medium. 



BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 presents a schematic related to a continuous-time implementation of the invention. 
It shows the key variables and how they combine to produce the value of the derivative of 
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the underlying variable x e . 

FIG. 2 is a flow diagram showing the steps of the method in the continuous-time case. 

FIG. 3 is a schematic diagram illustrating the projection of V + dV onto the space M. 

FIG. 4 shows a single step of a binomial model of variable A with probabilities pa and I — Pa- 

FIG. 5 shows a single step of a binomial model of a variable B when converted to risk-neutral 
form with risk- neutral probabilities of q& and 1 — q&. 

FIG. 6 presents a schematic of the method for the discrete-time case. It shows the key 
variables and how they combine to produce the value of the derivative of the underling 
variable B. 

FIG. 7 is a flow diagram showing the steps of the method in the discrete-time case. 



DETAILED DESCRIPTION: OVERVIEW 

Extended Black-Scholes equation 

Several embodiments of the invention are based on solving an extended Black-Scholes equa- 
tion. This novel equation differs from the standard Black-Scholes equation in that it involves 
a new coefficient in one of its terms. This coefficient allows the applications of the equation 
to be extended to include situations that the standard Black-Scholes equation fails to cover. 

A schematic depiction illustrating aspects of an embodiment of the invention is presented 
in FIG. 1. The figure shows an underlying variable x e , which is specified at 100. The 
variable x e defines a new asset payoff through a function F. Also shown in the figure is a 
market representative x m , determined at 110. The market representative x m is extracted 
from marketed securities, which are determined at 120. If this market representative is 
most correlated with x e , it will depend on x e ] and this is indicated in the figure by the 
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dashed line from x e to the market representation choice process, 110. The properties of 
this representative and x e define the particular coefficients of a partial differential equation, 
130. The function F defines the terminal boundary condition at 140. The solution of the 
equation then yields the value of the new asset at 150. 

Although the new method may at first appear arbitrary since it is based on projection, 
it has a compelling justification. The Black-Scholes approach uses replication to render the 
new payoff redundant, in the sense that the payoff is already embedded in the securities 
market. In a similar way, the new approach approximates the payoff and renders the payoff 
irrelevant in the sense that no risk-averse investor will want to own the new payoff because 
it is inferior to things already embedded in the securities market. 

Optimal hedge 

Once the value function V(x e ,t) is found, it is possible to find the portfolio strategy that 
best approximates the new asset, and can be used as a hedging strategy. 

When the hedging portfolio has value H(x e ,t), an amount (H(x e ,t) - <j>) is invested in 
the risk free asset and an amount (f> is invested in the most-correlated marketed asset where 

</>{x e ,t) = V x (x e ,t)x € p ec . (3) 

The minimum variance of the error V(x e , T) - H{x e ,T) can be found by solving another 
partial differential equation, namely, 

S t (x e ,t) + S Xe (x e) t)fi e x e + \S XeXe (x e ,t)olxl + e 2r(T ^[\4 e (^ e ,t)a e x e ] 2 (l - p 2 ec ) = 0, (4) 

with boundary condition S(x e) T) = 0. The value S(x e , 0) is the variance of the replication 
error at T, as seen at t = 0. 



LUB-101/US 



8 



Risk-Neutral Form 



An embodiment of this invention is an implementation as a practical risk-neutral valuation 
method. In this version the drift of the underlying variable x e is changed from /z e to [jj, e - 

Pernio™ ~ r)}. 

The general formula for value is then 

V(x e ,t) = e- r ^E[F(x e (T))}, (5) 

where E denotes expectation with respect to the risk-neutral model. Hence, in this version, 
the value is simply the discounted expected value of the final payoff, with expectation taken 
with respect to the risk-neutral model. 

The advantage of this form of the valuation is that value can be estimated by simulation. 
Basically, thousands of simulation runs of the underlying variable are made and the average 
discounted value is used as the expected discounted value. Such simulation methods are 
standard practice. A new feature of this embodiment is specification of the risk-neutral 
equation. 

Discrete-Time Version 

The risk-neutral form can be adapted to a discrete-time framework where the dynamics of 
random variables are represented by lattices. In the basic model, the underlying variable 
B is assumed at each step to either go up (with probability p B ) or down (with probability 
1 - p B ). A payoff that is function of this underlying variable can be evaluated by using the 
risk-neutral probabilities, denoted qg and 1 — q B . These, in turn, are found by using a market 
representative as in the continuous-time case. If this market representative is denoted A with 
price (value) va, the risk- neutral probability q B is 
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q B = Pb-cov[A,1(U b )][A-v a R}/o 2 a (6) 

= PB - Pl(U B ) t A& - V A R ] W 

where 1(U B ) is a payoff of 1 if B is up, and 0 if B is down. The expression "cov" stands for 
covariance, and Pi(u b ),a is defined as cov(l(tf B ), j4)/variance(A). 

A derivative of B is a payoff with a value of G u if up occurs and payoff G d if down occurs. 
The current price of such a (future) payoff is 

V B = ±[q B G u + (l-q B )G d }. (8) 

This is the basic method for the discrete-time case. The unique feature is the formula 
(6) for the risk-neutral probability. 

Recursive Solution and Hedging 

If the payoff occurs at the end of some finite number of steps (at the end of a lattice structure 
for 5), the value function can be found recursively. The recursive solution is 

= ^[q B V k (sU) + (1 " lB)V k (sU)l (9) 
where denotes the upper successor state to s*_i and sf^ denotes the lower successor 
state to 5 fe _i. The process is started with the terminal boundary condition specifying the 
payoff of the derivative G. 

An optimal hedge H can be found. If at the beginning of time period k the hedge 
portfolio's value is F fc -i, then H k -i - ^V k -i is invested in the risk free asset and -yV k ^ x is 
invested in the most-correlated asset A, where 

7 = cov(V t /V fc _ 1> i4)M. (10) 
The variance associated with optimal hedging can also be found by a recursion, as discussed 
later. 
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A Complete Methodology 

Embodiments of the invention provide complete methodologies for treating payoffs that 
are derivative of non-traded variables. It has a strong justification and a straightforward 
representation. For many problems, the new method can be put into practice by augmenting 
standard methods of financial computation. The requisite equations or lattice structures are 
similar to those derived from the Black-Scholes methodology and hence the computational 
methods are similar. These are methods for solving partial differential equations, methods 
for evaluating lattices, and simulation methods. A significant additional step of the new 
method is determination of the market representative x mj and estimation of the constants 
fj, e , fi m (the drifts), the variance of the market representative and the value of /? cm . 

Applications 

Embodiments of this invention have applications in numerous financial areas. For example, 
the methods may serve four functions. The first is pricing, as in determining the fair value of 
an option on a non-marketed variable. The second is design, as in the design of optimal con- 
tracts, project plans, business arrangements, and various agreements. The third is hedging, 
to minimize the risk associated with a non-traded asset by use of offsetting market partic- 
ipation. The fourth function is that of risk assessment, determining the residual risk after 
optimal hedging. In all cases, a financial payoff is tied to the performance of a variable that 
moves randomly but which is not a marketed financial security. The variable that defines 
the payoff is termed the underlying variable. 

The key idea is that when the underlying variable of an asset is not marketed, and hence 
the standard Black-Scholes equation does not apply, a market representative can be used 
instead. 

One large area of application is to the pricing and hedging of "off-exchange" derivatives. 
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For example, it may be desirable to hedge grapefruit production. Here the underlying vari- 
able is grapefruit price but there may be no financial instrument (such as a grapefruit futures 
contract) that is directly related to this underlying variable. Hence grapefruit production 
is not a derivative security and the Black-Scholes theory is not applicable. Nevertheless, a 
combination of orange juice futures contracts may be most-correlated with grapefruit prices 
and thus can serve as the market representative. This representative can be used (according 
to the valuation equation) to find value and also serve as the basis for an optimal (but imper- 
fect) hedging strategy. Similar situations arise with off-exchange trading of energy contracts, 
agricultural products, metals, foreign exchange, and many others. 

One important example is that of real estate projects. Typically, there is no underlying 
security for a real estate project, and thus the Black-Scholes equation does not apply. How- 
ever, the payoff of such a project is likely correlated with interest rates, real estate investment 
trusts, and regional economic variables. A combination of these can serve as the appropriate 
market representative and hedging strategy. 

Another general area of application is the determination of long-term futures prices. 
Typically, futures contracts extend several months, or at best a few years. Yet, many large 
projects span several years and correct pricing and hedging is important. The methods of 
this invention can be used to obtain such prices and hedges. 

Bonds are often issued that have unique risk, not completely correlated to marketed 
securities. For example, a corporate bond's price may be only loosely correlated to interest 
rates and to the stock price of the associated firm and to the stock prices of firms in the 
same or related industry. Hence the bonds are not perfect derivatives, but the methods of 
this invention can be used to evaluate them. Large physical projects, such as oil rigs, dams, 
various other infrastructure projects and private or public works have cash flows that are 
to some extent associated with risks that are not derivatives of market securities. Projects 
in foreign countries have associated country risk, which can only partially be hedged with 
bonds and insurance. 
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Contracting is another large area of application. For example, an electronics firm may 
plan to purchase a great deal of DRAM during the year. The firm can design contracts with 
options and price caps that provide reduced risk at favorable prices. The design of such a 
contract can be deduced from the market representative and the methods of this document. 
Contracts for movie and publishing rights can also be designed this way. 

The purchase of intellectual property rights through licenses is a fertile area of applica- 
tion. For example, a pharmaceutical firm may wish to license a biochemistry patent from a 
university, and both the design and the pricing of such a license could be carried out with 
the methods. In this case, the market representative may be a combination of the security 
prices of other pharmaceutical companies, HMOs, and other insurance companies. 

Insurance and various guarantees generally embody a great deal of risk. Such risk can 
be priced and hedged by the insurance company or guarantor by these methods. The risk of 
fire in an area may, for example, be correlated with weather derivatives and energy prices. 

Mortgages contain interest rate risk, but also prepayment risk. The prepayment risk 
may be correlated (imperfectly) to stock market indices, real estate prices, as well as interest 
rates. 

Real options is a general term which includes options related to many of the examples 
listed above, but especially to business projects. The methods of this document can improve 
both the evaluation and design of such projects by using a theoretically justified valuation 
formula rather than an ad hoc use of the Black-Scholes equation. 

Another important area of application is that of overall risk assessment. The residual 
risk of each individual project or derivative can be determined by the variance given as the 
solution to a partial differential equation that is used in embodiments of the invention. The 
risk associated with a collection of such projects can be found by solving a similar equation 
that depends on the value function for each of the projects separately. 
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DETAILED DESCRIPTION: 
CONTINUOUS-TIME CASE 



This section presents example application areas and gives a detailed description in both the 
continuous-time and discrete-time frameworks. 

Outline of Steps: Continuous-time framework 

An outline of the steps used in one embodiment of the invention are indicated in FIG. 2. 

1. Set Up. In this step, 200, we identify an asset of interest which gives a future payoff 
(or payoffs) according to the value of some variable that moves randomly with time. 
This variable is termed the underlying variable. If this variable is itself traded (and is 
governed by geometric Brownian motion (GBM)), then the asset of interest is a pure 
derivative and the standard Black-Scholes equation can be used. If the underlying 
variable is not traded, then the extended Black-Scholes equation is used. 

For example, an asset might be the future harvest of a grapefruit farm. The underlying 
variable determining payoff is grapefruit price which varies randomly from the current 
time until harvest. This underlying variable cannot be traded. It is true, of course, 
that grapefruit can be bought and sold at any time, but grapefruit cannot be held as a 
security since there are storage costs and grapefruit is perishable. The grapefruit price 
is denoted x e and is assumed to follow a GBM process of the form 

dx e = ti e x e dt + a e x e dz e (11) 

where z e is a standardized Wiener process. The process starts at time 0. The asset to 
be priced has a payoff at time T of F(x c , T). 

There are n (perhaps thousands of) securities that can be traded. Each of these follows 
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a similar processes, namely, 

dxi = faXidt + GiXidzi (12) 

where x { is the value of the i-th asset, for 1 = 1,2, ... ,n. The standardized Wiener 
processes Z{ (of zero mean and unit variance) are correlated, with cov(dZi, dzj) = Pijdt. 
Among these n assets there is an asset that is risk free at each instant. By forming 
an appropriate combination of the marketed assets, it may be assumed without loss of 
generality that the risk free asset is x n . Then 

dx n = rx n dt. 

In the standard Black-Scholes method, only the underlying security, the risk free rate 
r, and the payoff function is specified. No consideration is given to other securities, 
because the underlying security is all that is needed. In embodiments of this invention, 
other securities are considered in order to select a market representative, as explained 
in the next step. 

2. Determination of a Market Representative. The extended Black-Scholes equation uses 
a market representative x m , which is a combination of marketed assets. In this step, 
210, this market representative is determined. Generally, this market representative 
can be chosen in one of several ways. 

(a) A most-correlated asset is a marketed asset x c (or a combination of marketed 
assets) whose instantaneous return is most correlated (or close to most correlated) 
with dx e . In theory, it is obtained by solving the equations 

n-l 

£ VijVj = ^ie (13) 

for each z = 1, 2, . . . n — 1. The resulting a;'s are then normalized to sum to 1. 
The most-correlated asset is the asset x c with instantaneous return 

dx n ~ 1 dr 

^l = J2a^, (14) 

X c j_i X{ 
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x c can be used as the market representative x m . 

In practice, a most-correlated asset is obtained by considering only those securities 
that are obviously closely related to the underlying variable x e . 
In the case of grapefruit, for example, a most-correlated asset might be a futures 
contract on orange juice. (There are no grapefruit futures.) Such contracts are 
marketed securities, and are likely to be highly correlated with grapefruit prices. 
As another example, in seeking a most-correlated asset for DRAM (computer 
memory chips), one might consider the stock of several companies that produce 
DRAM and then find the combination of these stocks that is most correlated with 
DRAM prices. 

Ideally, this step entails identification of suitable candidates, evaluation of the co- 
variances of these candidates among themselves and with the underlying variable, 
and optimization of the combination of these candidates to provide the highest 
correlation. These steps are not difficult provided one has historical data on the 
security prices and the underlying variable. 

However, it should be understood that correlations, variances, and drift rates are 
measured approximately, and hence in practice exact parameters are not known, 
even for the standard Black-Scholes equation. Likewise, it is impossible to prove 
that a given asset is, in fact, the asset that is 'most' correlated with the under- 
lying variable. As in any method driven by estimates from limited data, one 
attempts to balance accuracy with data gathering and processing effort. Thus, 
in the context of the present description, a 'most correlated' asset or variable is 
defined to include not only the maximally correlated asset or variable, but also 
approximations thereto. 

Although in practice a 'most-correlated' asset is not known with certainty, this 
embodiment is distinguished from the following one by the fact that the search for 
the strongest correlation is restricted to obvious candidates within one or a few 
related market sectors or industries. For example, the steel industry likely would 
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not be considered as a candidate for an asset of grapefruit production. This 
contrasts with the Markowitz portfolio (generally approximated by the market 
portfolio) which includes all marketed securities. 

(b) The Markowitz or market portfolio. This alternative embodiment produces the 
same price as the method described above. This method is typically less conve- 
nient than the first and does not lead to an optimal hedging (or, equivalently, 
replicating) portfolio. In theory, this alternative is determined by solving the 
system of equations 

n-l 

These a^'s are then normalized to sum to 1 and used to define the portfolio with 
instantaneous return 

dx M ^ dxi _ 

= 2L, = i^Mdf + <JmVZM< 

Then xm can be used as the market representative x m . 

As an approximation to the Markowitz portfolio, the market portfolio consisting 
of the capitalized weighted average of the marketed assets may be used. For 
example, a diversified mutual fund can be used, or the S&P 500 index. Each of 
these is nevertheless only an approximation. It would be essentially impossible to 
determine the true Markowitz portfolio because of the inherent error associated 
with determining the required parameters. 

In contrast to the previous embodiment, the steel sector, and indeed all sectors 
of the economy, would be represented in this method. The advantage of this 
method is that the Markowitz (or market) portfolio, once found, can be used for 
derivatives of any underlying variable. 

(c) Local index. Suppose there are two non-marketed variables x ei and x e2 that deter- 
mine various derivative assets. Let x c i and x c2 be corresponding most-correlated 
risky marketed assets. Next let x m be the Markowitz combination of x cl and x c2 . 
Then the single asset x m can be used as a market asset for all derivatives of x e \ 
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and x e 2- This procedure can be directly extended to any number of underlying 
variables. 

To complete the set up, the payoff function F(x e {T)) ) the parameters \i e and a e of the 
underlying variable, the parameters /z m , /J em , and the interest rate r are specified. As 
stated earlier, in practice it is not necessary to specify the parameters of every marketed 
security. A market representative can be determined by consideration of those assets 
closely related to the underlying variable or, alternatively, to a broad index of the 
market. 

3. Define the Differential Equation. In step 220 various terms are specified for the ex- 
tended Black-Scholes equation: 

rV(x e ,t) = V t {x e ,t) + V Xe (x e) t)x e [fi e - f3 em {fi m - r)] + lV XeXe (x ey t)x 2 e al (15) 

where 

Pern = E[dz e dz m ] 
Pern = Pem&e/Gm- 

The subscript m denotes the market representative. The boundary conditions are 

V(x e) T) = F(x e (T)) 
V{0,t) = e- r(T -^F(0). 

Equation (15) plays a central role in several embodiments of the invention. It dif- 
fers from the standard Black-Scholes equation in the coefficient [/z c - f3 e m{pm — 
This coefficient replaces r in the standard Black-Scholes equation as a coefficient of 
K e (z e > The new coefficient uses properties of the market representative x m . Equa- 
tion (15) reduces to the Black-Scholes equation if the new asset is a pure derivative so 
that x m = x e is the most-correlated market representative. 
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4. Solve the Differential Equation. In step 230 the differential equation is solved. This 
can be done analytically in some cases. Most often it will be solved by standard 
numerical techniques. Once the parameters are known, solution of the equation is no 
more difficult than solution of the standard Black-Scholes equation. It is only one 
coefficient that is different. 

The value of V(x e (Q),Q) is the proper price of the derivative. The procedure can be 
terminated here, with the proper price. Additional information is obtained in the next 
two steps. 

5. Find the Optimal Replicating Portfolio. Optionally, step 240 calculates the best repli- 
cating portfolio. One may replicate the market portion of the risk in the asset by hold- 
ing an appropriate portfolio of market securities with initial value equal to V(x e (Q),0). 
This portfolio is a combination of the most-correlated asset and the risk free asset. 
The amount to be held in the most-correlated asset as a function of x e and t is 

<t> = V Xe (x e ,t)x e (3 ec 

where in this case, we use the subscript c to emphasize that the market representative is 
taken to be the most-correlated asset. Holding the negative of this replicating portfolio 
produces an optimal hedge for the derivative. This procedure is essentially identical 
to that used in standard derivatives of the Black-Scholes type. The difference is that 
in the Black-Scholes case the hedge is perfect, while in embodiments of this invention 
the hedge may not be perfect. 

For example, in the case of grapefruit production, the optimal hedge would be a port- 
folio consisting of the risk free asset and a position in orange juice futures. The propor- 
tions of the two components of the portfolio would be adjusted frequently according 
to the formula for </>. If the hedge is in place, the hedge portfolio yields cash immedi- 
ately equal to V(x e (0), 0) and the final net payoff (consisting of sale of the grapefruit 
production and settling of the hedge contract) is on average zero, but it has a residual 
variance. 
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Ideally, the portfolio adjustment would be continuous, but in practice periodic adjust- 
ments (daily, weekly, or monthly) are sufficiently accurate. These types of adjustments 
are common in pure derivatives methodology as well, where the ideal hedge is perfect. 

6. Find the Error Variance. Optionally, step 250 calculates the tracking error. The 
variance of the part of the payoff that cannot be hedged is found by solving another 
partial differential equation which uses the solution to (15) as an input. The equation 
is 

S t (x e ,t) + S Xef i e x e + \S XeXe alx\ + e 2 < T ^[V Xe a e x e ] 2 (l - p 2 ec ) = 0 (16) 

with boundary condition S(x €7 T) = 0. The value of S(x e (0),0) is the variance of the 
error at time T as seen at time 0. This may be extended to computation of the variance 
associated with combinations of derivatives. 

For example, in the case of grapefruit, the value of S(x e (0),0) gives the variance of 
the exposure at time T with an optimal orange juice futures hedge. That is, with the 
hedge, V(x e (0) ) 0) is obtained at t = 0 and the net amount (of the combined grapefruit 
production and hedge) attained at t = T has expected value 0 but variance equal to 
S(x e (0),0). 

This calculation of error variance has no analogy in the standard Black-Scholes case 
because in that case the error is always zero (in theory). Hence, this method for finding 
the hedge variance is new and important. It is not difficult to solve the appropriate 
partial differential equation, for indeed, it is similar in structure to the equation for 
value. 

Extensions 

There are several important extensions of the method. For example, 

1. Inclusion of cash flows that occur with time, with incremental cash at time t being of 
the form h(x ei t)dt. 
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2. Varying parameters. The parameters of the model can vary with t, x e , and the x^s. 
In general this leads to a partial differential equation of higher order. 

3. Parameters, such as the risk free rate, may be governed by stochastic processes. 

4. Additional variables. Additional non-marketed variables can be introduced. These 
serve as "state variables" for the system. For example, a state variable might be total 
industry productive capacity, or estimates on the probability that certain legislation 
will be passed. 

5. Estimation variables. A suitable non-marketed but observed variable may be the esti- 
mate of an unobserved variable that serves as the underlying variable for the payoff. 
If the estimate converges to the actual value at the time of payoff, the estimate may 
be used at all points instead of the original variable. For example, the best estimate 
for yearly revenue may converge to the actual figure as the year ends. 

6. Non-market random components. If the final payoff is a function of the form F* (x e , y , T) 
where y is random and independent of x e and the market, then we define F(x e ,t) = 
E y [F*(x e ,y,T)} where E y denotes expectation with respect to y. The value of the 
payoff can then be found with this F and it will be the proper price. In this case, 
the boundary condition for auxiliary differential equation for variance is S(x e ,T) = 
variance y [F* (x e , y , T)] . 

7. Alternative Processes. The method of relating the non-market variable to the market 
in order to obtain the correct price can be extended to alternative processes, including 
jump processes, and moving average processes. 

8. Market Cash Flows. In some cases the cash flows or payoffs may depend on marketed 
variables as well as non-marketed variables. In such a case the value function V will 
depend on both x e and the marketed variables. The most-correlated asset for the 
marketed variables will be those variables themselves, while the market representative 
for the non-marketed variable is found as in the basic case. 
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9. Path Dependent Cash Flows. These are cash flows such as max t [0,z c (£)], 0 < t < T 
that depend on the actual path taken by x e rather than on its instantaneous or final 
value. The risk-neutral process can be used to evaluate such situations by taking the 
risk-neutral discounted expected value of all cash flows. 

Brief Derivation of the Equations 

We refer to step #1 (Set up) to define the problem. Hence, we have an underlying variable 
x e governed by geometric Brownian motion (GBM) over 0 < t < T as 

dx e = fj, e x e dt + a e x e dz e) (17) 

and likewise there are n marketed assets that also are governed by GBM. There is a special 
asset of interest whose payoff at time T is F(x e (T)). 

We propose an (initially unknown) price function V(x e ,t) on the time interval [0,T] with 
terminal value V(x e ,T) = F(x e ). It will be defined by instantaneous discounted projection. 

At a fixed time t we define the space M as the linear space generated by all instanta- 
neous marketed returns (dx^/xi. That is, M is made up of all linear combinations of those 
instantaneous returns. Symbolically, the value of our V(x e , t) will satisfy 

V{x e ,t) = P{V{x e ,t)+dV{x e ,t)\M}. (18) 

which denotes the discounted projection P onto the space M of instantaneous marketed 
returns. Specifically, V(x e ,t) is found by: first calculating the payoff m* at time t + dt in M 
that is closest to V(x e , t) + dV (x e , t) in the sense of minimizing the expected square of their 
difference; and second, computing the price at time t of m* (which is defined by linearity 
since it is a combination of marketed returns). Briefly, we say that the value at t of the 
payoff at t + dt is given by discounted instantaneous projection. This process of projection 
onto M is illustrated schematically in FIG. 3. 
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Since payoff V(x e ,t) at t + dt is certain, its price at t is the discounted value (1 - 
rdt)V(x e ,t). Thus (18) becomes 

rV{x„t)dt = P{dV(x e ,t)\M}. (19) 

We shall only keep terms that are first-order in dt. Hence, in particular, P{dt} = dt. 
Substituting the Ito formula 

dV{x e ,t) = [VUx,,«) + V«.(x ei t)^ (20) 

into (19) and keeping only first-order terms in dt, (19) becomes 

rV {x e , t)dt = [V t (x e: t) + V Xe (x e) t)fi e x e + lV XeXe a 2 e x 2 e }dt + V Xe o e x e P{dz e \M). (21) 

The projection of dz e is a linear combination of rdt and the differential of a marketed 
asset most correlated with dx e . A marketed asset x c with instantaneous return dx c /x c most 
correlated with dx e /x e satisfies drr c = /vc c d£ + a c x c dz c where \i c and a c are each formed as 
(identical) linear combinations of the ^'s and <Vs, respectively. z c is a standardized Wiener 
process. 

Standard methods show that the projection of dz e onto M is p ec dz c . In addition, from 

x c = P{x c + ii c x c dt + x c a c dz c } 

it follows that 



P{ dZc } = (r Mc) dt. (22) 



Substituting this in (21) produces 



rV{x e) t) dt = [Vt(x c ,t) + V Xe [x ei t)fi e x e + |K eXe (x € ,t)a^]dt + V Xe P{p ec a e x e dz c } 

= [V t (x„t) + V Xe (x e ,t)x e [ll e - /? ec ( Mc - r)] + iKeXe^e^)^]^, 



where we define 
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Canceling the dt we have the final result 

rV(x e , t) = Vt(x e , t) + V Xe {x e , t)x e [fie - Pec&c ' r)] + \V XeXe {x e) t)o\x\. (24) 

This is the basic extended Black-Scholes equation. 

Once equation (24) is solved with the boundary conditions V(x e ,t) ~ F(x e ) ) V(0,t) — 
e" r(T "^F(0), the value V(x e (0),0) is the value at time 0 of the derivative. 

Market Representatives 

Details of how the two basic alternative market representatives, a most-correlated asset and 
the Markowitz portfolio, are computed are given here. In practice, these detailed calcula- 
tions are greatly simplified (approximated) by use of intuition or equilibrium arguments, as 
mentioned below. 



1. Most-Correlated Asset. 

Arrange the n marketed assets so the first n — 1 are risky and the n-th is risk free. The 
return of a most-correlated asset is of the form 

dr n ~ l dr- 

H£ = £oA (25) 
where the a^s do not necessarily sum to 1. 

To maximize the correlation of this return with dx e , the c^'s solve 

max E 



/ dx c dx e \ 
\x c x e J 



( dx c 

subject to var I — - ) < 1. 



x c 

Introducing a (positive) Lagrange multiplier for the constraint, and representing dx c /x c 
in terms of the c^'s, this is equivalent to solving 

n— 1 n—1 



max 



ot,i<j ie - x J2 a i&ij a j 



(26) 
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Since any solution can be scaled by an arbitrary positive constant and still preserve 
the property of being most correlated, without loss of generality one may set A = 1. 
Then (26) becomes 



which is easily solved. 

The current price (at t) of dx c /x c is c^. To make the price equal to 1, the c^'s are 
scaled so that they do sum to 1. Then the correct x c is given by (25) with the scaled 



In practice, a market representative may be chosen as an asset that is obviously closely 
related to the variable x e . For example, an orange juice future is closely related to 
grapefruit price. It is probably not necessary to consider oil and metal futures or 
the thousands of marketed stocks to find a closely related market representative. A 
refined approach is to consider a family of agricultural futures and, based on historical 
returns, compute the most-correlated portfolio from this family. The selection of a 
closely related market representative is, in practice, partially art and partially data 
collection and computation. 

2. Markowitz portfolio. The Markowitz portfolio xu has the advantage that it can be 
used to price a derivative of any process x e , but its disadvantages are that is difficult 
to use, less intuitive, may not exist, and does not lead to an optimal replication. (The 
use of this portfolio in the single-period case is described in Luenberger [10].) The 
Markowitz portfolio x M is the portfolio of purely risky assets (the first n — 1 assets) 
that has price 1 and maximizes 

If it exists, the set of a/s that achieves it satisfies 



^2 VijUj = °ei 
3=1 



(27) 



CXi'S. 



n-1 



J2 OijCij = c{m - r) 



(28) 
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for the constant c that makes the a^s sum to 1. This portfolio is an alternative for 
purposes of valuation only, not for determining optimal replication or hedging. In 
practice, a broad market index, such as the S&P 500, is used as a proxy for the 
Markowitz portfolio. 

3. Local index. Suppose x ei and x e2 are non-marketed underlying variables and x ci and 
x c2 are the respective most-correlated marketed assets. The Markowitz combination 
of x c i and x c2 is the portfolio x m = a x x c i = a 2 x c2 where the weights ot\ and a 2 are 
determined as in the above discussion of the Markowitz portfolio under the assumption 
that x c i,x c2 and r are the only marketed assets. This x m is a market representative 
for any derivative of rc el and x e2 . It serves as a local index, and may in practice be the 
index of an industry to which x ei and x e2 are related. This local index can be used 
to value a variety of derivatives. It cannot be used as the basis for optimal hedging, 
but it can be used to form approximate hedges. The method extends to any number 
of underlying variables. 



Important Special Case. 

If the payoff is a put or call option on the variable x e , a closed- form expression for the 
solution to (15) applies, which extends the special case solved by Black and Scholes. If the 
strike price of the call option is K, the payoff is max[x e — K, 0]. The value of the call option 
is 

V(x e , t) = xJ^^Nid,) - Ke" r ^N(d 2 ) (29) 



where N(d) denotes the value of the standard normal distribution at amount d, and 

+ {u> + \c 



* . M^^W-t) (30) 



d 2 = di - <r e VT - 1 (31) 

W = H e - PemifJ-m-r), (32) 

where p, m is the drift rate of the market representative x m and p em is the beta of x e and x m . 
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The value of a put with strike price K can be found from an extended put-call parity 
formula 

C-P + Ke~ rT = x e e^ T . (33) 
Here C is the price of the call, and P the price of the put. 

Optimal Replication 

It is possible to optimally replicate the final payoff by trading in the marketed assets. The 
replication is imperfect, but the error has zero expected value and is uncor related with 
all marketed assets. This is, therefore, the best that can be done within a mean-variance 
framework. The trading strategy for replication can be derived from the value function. 

As in (20), the value function follows the Ito process 

dV(x e ,t) = [V t {x e ,t) + V Xe {x e) t)fi e x e + ±V XeXe ^ (34) 

One may write 

dz e = p ec dz c + yjl- p 2 ec dz P) (35) 

where dz p is a standardized Wiener process uncorrelated with all marketed assets and 
where p ec = o ec /{a e (r c ) = Pec&c/Ve- With this substitution and substitution of the term 
\V XeXe (x e ,t)alx 2 e from the extended Black-Scholes equation (15) it is possible to transform 
(34) to 

dV(x e) t) = [(V(x e , t) - <t>)r + <t>n c ]dt + </>a c dz c + ddz p , (36) 

where 

<j>{x e ,t) = V x {x e ,t)x e p ec (37) 
6(x e ,t) = V Xe {x e ,t)a e x ey /l-pl. (38) 

Motivated by (36), the process H is defined by 

dH{x e , t) = [(H(x e , t) - <j>)r + <pii c ]dt + <t>a c dz ci (39) 



LUB-101/US 



27 



with initial condition tf(x e ,0) = V(x e ,0). 

From (39) it is clear that dH (x e , t) is in M at (x CJ £) and that, overall, #(x e , t) is generated 
by linear combinations of marketed returns. From 

d(V - H) = r(V - # )d* + 5dz p . (40) 

and the initial condition V - H = 0, it is clear that at time 0 there holds E[V(x e> t) - 
H(x e ,t)} = 0 for ail £, 0 < t < T. Furthermore, since V - H is a linear integral of d^'s, it 
is orthogonal (that is, uncorrelated) with all marketed returns. Thus the random variable 
H(x ei T) as seen at time 0 is the projection of V(x e ,T) onto the space of marketed assets 
over 0 < t < T. 

It follows that the optimal replication is simply H, governed by (39). This replication 
requires only initial cash of V(x e ,0) and no other infusions or withdrawals. At each instant 
H — (j) is invested in the risk free asset and <j> in the most-correlated asset (for a total of H) . 
Holding — H serves to optimally hedge the derivative. 

Projection Error 

The variance of the residual replication error V - H at T can be found as a solution to a 
partial differential equation adjunct to the general pricing equation. 

Define the difference variable D = V ~ H , and let U = e r ^ T ^D. Then from (40) 

= {-re r ^D + e T{ ?-%D}dt + e r ^6(x e , t)dz p 

- 5*{x e ,t)dz p , (41) 

where 6*(x ei t) = e r ^ T '^S(x ei t). 

Let S(x e) t) be the variance of U(T) when in (41) U is initiated with U = 0 at the point 
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(x e ,t). This variance is 



S(x e ,t) = E J 6*(x e (s),s) 2 ds , 



(42) 



where E denotes expectation at t. 



We write (42) as 




- 8*(x e , t) 2 dt + E [S(x e , t) + dS{x e , t)} . 



(43) 



Hence. 



E{dS{x e ,t)} + 6*(x e ,t) 



\ 2 dt = 0. 



(44) 



Expanding (44) using Ito's lemma, we find 



St(x e , t) + S Xef , e x e + \S x ^a\x\ + e 2r ^[V Xe (x € , t)a e x e ] 2 (l - p 2 ec ) = 0, 



(45) 



with boundary condition S(x e ,T) = 0. The value S(x e ,0) is the variance of the replication 
error at T, as seen at t = 0. 

Universality 

As a by-product of the optimal replication equations, we obtain a compelling justification 
for methods used in embodiments of this invention. When the new asset is priced according 
to teachings of this invention, the optimal replication H has the same expected value as the 
asset with payoff F{x e (T)). In addition the optimal replication has lower variance than the 
asset since V(x ei T) = H(x e ,T) + error, where the error is uncorrelated with all marketed 
assets. This means that every risk-averse investor will prefer the replication over the new 
asset. This property is termed universality, to reflect that it is true for everyone. Another 
way to state this property is that if the new asset is priced according to methods of this 
invention, every risk-averse investor will elect to include the asset only at the zero level 
(preferring neither to be long or short the asset). In this sense, the price renders the asset 
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irrelevant in the market, for it is dominated by the best replicating asset which is already 
available. 

Risk-Neutral Process. 

These methods can be transformed to formulas based on a risk-neutral processes. The novel 
feature is that the risk-neutral process is defined in terms of a market representative x m 
defined earlier. The appropriate risk-neutral process is 

dx e = ux e dt + a e x e dz e (46) 

where 

LJ = // e - (3em{V>m ~ r). (47) 

In these terms the proper price of the payoff F(x e (T)) is 

p = e- rT E[F(x e (T))] (48) 

where E denotes expectation at time 0 with respect to the risk-neutral process. One way 
to solve such a problem is by simulation. Many paths of x e are generated according to the 
risk-neutral process and the resulting terminal payoffs are averaged to obtain an estimate 
of the risk-neutral expected value. The discounted value of this estimate converges to the 
correct price as the number of simulation runs increases. 

A significant advantage of the risk-neutral process is that it can be used to value path- 
dependent cash flows; where the payoff depends on the particular path taken by x e as well 
as its final value. For example, the Asian option with payoff max [max 0 <t<T £ e (£) — K, 0] can 
be evaluated using the risk- neutral process. The proper price is again 

p = e- rT E[F(x e (T)}, 

which can be evaluated most easily by simulation, but also by special grid methods. 
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DETAILED DESCRIPTION OF THE INVENTION: 
DISCRETE-TIME CASE. 



Approximation 

One way to work in a discrete-time framework is to directly discretize the model and the 
partial differential equations. For example, the new version of the process for x e becomes 

x e (k + 1) = (1 + fi e At)x e + a e x e e(k)\/ r Ki. (49) 

In this equation At is the length of the time step of the model and e(k) is the value of a 
standardized normal random variable uncorrelated with previous or future such variables. 
The market variables are transformed to discrete form in a similar way. A corresponding 
discrete-time version of the extended Black-Scholes equation can be developed. It is 

- 1)) - ± {E[V k{Xe{ k)} - cov(H(x eW ,, m W)(E[, mW ]-^) | (5Q) 

where R = e tAt and where x m is the market representative (either the Markowitz portfolio 
or a most-correlated market asset or a combination of these.) The expected values are taken 
at time k — 1. The quantity p m is the price at time k — 1 of x mi is the variance of x m (k) 
as seen at k — 1. The formula is valid for any discrete process that governs x e , although the 
market variables follows a fixed process as in (49) (but with i's instead of e as subscripts). 
The extensions applicable in the continuous-time framework are also applicable in discrete 
time. 

The discrete-time approximation (49) can be converted to risk-neutral form as 

x e (k + 1) = (1 + uAt)x e + a e x e e(/c)\/A*, (51) 
where u = fi e - /? em (// m - r). 
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Finite- State Models 



A single stochastic variable A is represented by the two-state model shown in FIG. 4. The 
U A outcome is considered an "up" move, and the D A outcome a "down" move. The model 
has three degrees of freedom (one probability and two node values) and hence the expected 
value and variance of the variable can be matched with one remaining degree of freedom. 
The model can be extended over several time periods as a binomial tree or (frequently) a 
binomial lattice. 

There are two standard approaches to parameter matching: additive and geometric, cor- 
responding to matching moments of the variable itself or matching moments of the logarithm 
of the payoff. Geometric matching is a natural choice for processes governed by geometric 
Brownian motion. Additive matching is frequently used for discrete-time models. A binomial 
model for one of these is easily converted to a corresponding binomial model for the other. 
For small At defining the time step of the model, the two approaches are nearly identical. 

Suppose that A and B are variables as above, and suppose G is a function of (a derivative 
of) B. G is defined by the two values G u and G d corresponding to whether U B or D B occurs. 
There is an important relation between the covariance of such a derivative with A and the 
covariance of B with A that holds when B is described by a binomial model. This is spelled 
out in the following result, easily proved by algebra: 

c ° v ^ G) =II^H' 4 - s) - (52) 

Suppose A is the market asset most correlated with B. The above result guarantees that 
in the case of binomial models, A is also most-correlated to all derivatives of B (since the 
two covariances are proportional). It follows that the three- variable model (with A, 5, and 
the risk free total return R) can be used to price all derivatives of B. 

Paralleling the development for the continuous-time case, A is a market representative 
and it may be taken to be a marketed asset most-correlated with the underlying B or, 
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alternatively, as the Markowitz (market) portfolio. 

Let G be a derivative of B, defined by its two values G u and G d . The projection price is 
(with an over-bar on a random variable denoting expected value) 

v G = ±[E(G)-cov(G,A)(Z-v a R)/<t 2 a ] 

= ^ [ P BG u + (l- PB )G d -[coy(l(U B ),A)G u + cov(l(D B ),A)G d }[A-v A R)/a A ] 
= | { [pb - Pi(u B) AA - vaR}} G u + [(1 - PB ) - Px { dM - v A R)] G d } , (53) 

where 1{U B ) and 1{D B ) denote payoffs of 1 if U B or D B , respectively, occurs; and where 

P HUb) , a = cov[1(U b ),A]/<t a . 

This can be written as 

V G - \ [q B G u + (1 - q B )G d ] (54) 

where 

Qb=Pb- (3i(u b ),a[A - v A R). (55) 

Expanding the beta term, we can calculate explicitly 

q B = Pb ~ Pb{1 ~ Pb)[A u - A d ][A - v a R]/<j% (56) 

where A u = E[A\U B ] A d = E[A\D B }. 

This is the basic method for the discrete-time case. The unique feature is the formula (55) 
or (56) for the risk-neutral probability of an up move. Once this probability is determined, 
it may be used in place of the true probability p A for the purpose of evaluating payoffs that 
depend on B. In other words, for purposes of evaluation the binomial model takes the form 
of FIG. 5. 

Recursive Solution 

The single-period structure can be extended to a multiple-period framework by piecing to- 
gether single periods, and this leads to a recursive solution. 
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Suppose the nodes of the lattice are numbered by the time index k and the state index 
s k which is the level of the node counting from the bottom. Briefly, we write the node as 
(fc, s k ). The value function of the derivative is a value at each node, and described as V k (s k ) 
at time point k and state s k . 

The recursive solution is 

= ^[q B V k (sU) + (1 - q B )V k {sU)], (57) 

where s k _ x denotes the upper successor state to s k _ x and s k _ x denotes the lower successor 
state to Sfc_i. The process is started with the terminal boundary condition specifying the 
payoff of the derivative G. If there are additional payoffs along the way, they are incorporated 
step by step in the usual manner. 

This method can be extended in the same ways as the continuous-time version. 

FIG. 6 is a schematic depiction of the method for the discrete-time case. This schematic 
has the same basic structure as that of FIG. 1 which depicts the schematic for the continuous- 
time case, with the exception that the variable V is found through a recursion process 600. 
The figure shows how the underlying variable B defines the new asset and how the market 
representative A is extracted from the market. The properties of these two variables define 
q B which defines the recursion, leading to the value of the new asset. 

Error Process 

Denote by {Vfc|M} the projection of V k on the market at time k — 1. This will be the 
projection of V k onto the space spanned by R and a market representative A most correlated 
with B, It is easily shown that 

{V k \M} = V k + cov(^, A)(A - A)/a\. 

Hence 

V k =V k + cov(14, A)(A - A)/a\ + e k , (58) 
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where e k is uncorrelated with the market. Using the pricing equation (50), one has 

RVk-i =V k + cov(I4, A)(R - A)la\. 

Eliminating V k) (58) can be written as 

V k = RV k . t - cov(V k , A)(R - A)ja\ + cov(V k , A) (A - A)ja\ + e k 

= [(l- 7 )JZ + 7i4]V fc _i+£ Jfc (59) 

where 

1 = cov(V k /V k _ u A)/a\. 

In the same manner as in the continuous- time case this shows how to select the best approx- 
imating (or replicating) portfolio H, with an amount H k -i — lV k -i in the risk free asset and 
7l4_i in the portfolio A (for a total of H k ^ x ). Thus H satisfies the recursion 

H k = (i/ fc _i - fV^R + + €fc . (60) 

Error Propagation 

The replication error D k = V k - satisfies 

D k = RD k + e k (61) 
where e fe is uncorrelated with the market and 

e k = V k -V k -cm{V k ,A){A-A)lo\. (62) 

It follows that 

var(e fc ) = var(V fc ) - cov(V t> Afja\. (63) 
In terms of the lattice parameters this becomes 



■{• 



var(e fc ) = {p A (l - p A ) - jfi^5tf% ^ [K ~ ^ (64) 
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If the dynamics of the variables A and B are stationary, the expression in brackets is constant, 
and hence the error variance is simply a constant times the square of the difference in the 
two successor V^'s. 

Let U k = R T ~ k D k . Then U k = t4_i + R T ' k e k . If S k is the variance of U T as seen starting 
at fc, then 

5 fc _!=:E[S fc ]-fi? 2 ( r - fc Var(6 fc ) (65) 

with terminal condition S T = 0. This can now be evaluated by a backward recursion in the 
lattice. 

Method Steps 

Following is a summary of the steps of the lattice method for the discrete-time case. These 
steps are depicted in FIG. 7 and parallel those for the continuous-time case illustrated in 
FIG. 2. 

1. Step 700 is the set up. Given a random payoff G that depends on a non-traded variable 
B that evolves randomly, formulate a binomial lattice model. 

2. Step 710 determines an appropriate market representative A as either (a) a market 
asset most correlated to B or (b) the Markowitz or market portfolio. 

3. Step 720 determines the risk-neutral probabilities by using the formula (56). 

4. Step 730 solves for the values of V on the lattice corresponding to the variable B using 
the recursion (57). Initiate with the terminal boundary condition V = G. The result 
at time 0 is the price. 

5. Optionally, step 740 determines the optimal replicating portfolio from H 0 = Vo, and 
thereafter at step k - 1 investing H k -i - jV k -i in the risk free asset and in the 
asset most correlated with B. 
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6. Optionally, step 700 determines the error variance by the backward recursion (65). 



IMPLEMENTATION 

The present invention may be implemented in various different ways on a computer. For 
the purposes of the present description, the term 'computer' is defined to include any elec- 
tronic digital information processor, including hand-held calculators, personal data assis- 
tants, pocket personal computers, laptop computers, desktop computers, and so on. The 
computer may be programmed by providing a computer readable digital storage medium 
containing instructions to execute various methods of the invention. The program or pro- 
grams are then executed on the computer. Information associated with the variable x m may 
be provided to the computer manually by a user, received from a data network, or retrieved 
from a storage medium into computer memory. The program then retrieves the information 
associated with the variable x m from the computer memory and performs calculations in 
accordance with various methods of the invention. The results of the calculations may be 
used to display or otherwise communicate information to the user, or they may be stored or 
transmitted digitally for further processing or later display. 

More specifically, the implementation of methods of the invention may take various dif- 
ferent forms. Following are representative examples of specific implementations. 

1. Numerical Solution. The straightforward way to implement the continuous-time method 
once the model is defined is by solution of the extended Black-Scholes equation (15). 
This can be carried out by standard numerical procedures, following the methods de- 
veloped for the ordinary Black-Scholes equation. The simple, one-dimensional version 
shown explicitly in (15) can be solved by a finite grid method: a large two-dimensional 
grid of points is defined with coordinates corresponding to x e and t and grid point values 
corresponding to the associated V values. In the simplest method, first-order deriva- 
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tives of V are formed as (normalized) differences in grid point values. The second-order 
derivative is the second-order difference of grid point values. Boundary conditions are 
typically imposed at t = T, and at lower and upper values of x e . In practice, if the 
problem is defined for all values of x e with 0 < x e < oo, an appropriate value for 
V(x e ,t) at a high upper boundary curve x e (t) is assigned as a boundary condition 
in addition to the terminal condition and the condition at x e = 0. It is possible to 
carry out the complete solution with a spreadsheet program such as Excel running on 
a general purpose computer. Experience has shown that the primary technical concern 
for simple methods such as this is that the time step between successive grid points be 
small in order to assure convergence. Small grid point widths in both x and t are used 
to get good accuracy. 

For higher-order equations associated with extensions of the basic method, it is conve- 
nient to use a professional software package or one written in a programming language. 
The value of a call option on a non-tradable variable, but with a stochastically varying 
interest rate has been solved by the inventor using the Excel Visual Basic package. 
The same techniques may be applied to solve for the variance according to the partial 
differential equation (45). In simple experiments, the same grid size as used for finding 
V worked for finding the non-hedgable variance. 

2. Explicit Solution. The explicit solution for call options (or for puts using the put-call 
parity formula) is easily solved by using either a table, or more likely, an approximate 
formula for values of the normal distribution. Again this can be carried out with Excel. 
It could be easily incorporated into a hand-held calculator. 

3. Discrete-time Version. The discrete time recursion is easily carried out by direct recur- 
sion. This method is relatively free from convergence issues, since both the model and 
the solution method are carried out with the same time steps (unlike the continuous- 
time case which is computed with a discrete process). 

4. Lattice Methods. The lattice method is easily implemented. In simple cases it may be 
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carried out with a spreadsheet program running on a computer. The method becomes 
more challenging to implement when the model parameters are not constant, and when 
higher-order lattices and trees are used. 

5. Risk-Neutral Computation. Beginning with a discrete model, the risk-neutral version 
of the processes can be implemented computationally by computing the risk-neutral 
expected value of the payoff function. This can be carried out with backward recursion, 
simulation, or by constructing a lattice that has the risk-neutral probabilities on its 
arcs and using backward recursion on that lattice. 

6. Simulation. Simulation is a powerful method for implementing methods of this inven- 
tion. Standard simulation packages can be used, such as a spreadsheet program like 
Excel, modest simulation packages such, as Crystal Ball, or more advanced statisti- 
cal packages. At a professional level, special software that accounts for the financial 
structure as well as advanced simulation concepts would be used. 

7. Optimization. In many situations, on-going decisions can be made that influence the 
value of the payoff. For example, an American-style option allows the owner to exercise 
the option at any time before expiration, and selection of this exercise policy is an 
important component of the analysis of the option. Similarly, in a business venture 
there are opportunities to expand, contract, delay, and so forth. These policies can be 
found with the methods of this document in conjunction with standard methods for 
policy optimization employed for ordinary derivative theory. 
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